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ABSTRACT 

This  paper  presents  the  results  of  using  the  constant  modulus  algorithm  (CMA)  to  recover  a  shaped 
offset  quadrature-phase  shift  keying  (SOQPSK)-TG  modulated  signal,  which  has  been  transmitted  over 
an  aeronautical  telemetry  channel  using  the  iNET  data  packet  structure.  The  iNET-packet  structure 
contains  known  data  bits  (the  preamble  and  asynchronous  marker  (ASM)  bits)  within  each  data  packet, 
which  can  be  used  to  determine  the  minimum  mean  square  error  (MMSE)  equalizer  and  is  then  used  as 
an  initializer  for  the  CMA  algorithm.  A  baseline  analysis  of  the  performance  of  the  MMSE-initialized 
CMA  equalizer  has  been  previously  done  using  average  bit  error  rates  (BER).  In  order  to  achieve  real¬ 
time  processing,  in  this  implementation  the  CMA  equalizer  weights  are  updated  using  a  new  data  packet 
for  each  adaptation  of  the  equalizer  during  which  the  bit  error  analysis  is  being  done.  Two 
implementations  of  the  block  processing  CMA  are  compared.  The  first  achieves  convergence  of  the 
equalizer  in  the  startup  phase  by  processing  a  single  fixed  packet  and  the  second  achieves  convergence 
by  processing  the  multiple  packets  of  received  data.  Perfonnance  evaluation  results  based  on  bit  error 
rates,  are  presented  for  these  two  methods  and  they  are  compared  to  a  previously  determined  baseline 
performance.  We  note  that  with  the  proper  choice  of  parameters  these  real-time  methods  can  achieve  the 
baseline  perfonnance. 


INTRODUCTION 

We  investigate  the  effectiveness  of  the  CMA  equalizer  in  recovering  an  unknown  data  bit-stream  that 
has  been  transmitted  over  an  aeronautical  telemetry  channel  with  shaped  offset  QPSK,  version  ‘TG’ 
(SOQPSK-TG),  using  the  iNET  data  packet  structure.  The  CMA  equalizer  has  been  widely  used  and  is 
the  most  popular  blind  adaptive  equalizer  in  use  today  because  of  its  relative  simplicity,  and  also  its 
good  global  convergence  properties  [1].  It  is  well  studied,  and  it  is  known  to  be  effective  for  signals  and 
constellations  that  possess  the  constant  modulus  property. 

Since  SOQPSK-TG  is  a  partial  response  continuous  phase  modulation  (CPM),  it  should  be  a  perfect 
candidate  for  this  equalizer,  but  previous  work  by  Law  [2]  had  noted  some  difficulties  with  the  use  of 
CMA  to  equalize  SOQPSK-TG  over  a  telemetry  channel.  Previous  work  in  [3]  showed  that  a  CMA 
equalizer  with  block  processing  is  an  effective  equalizer  for  recovering  SOQPSK-TG  modulated  signals 
in  the  the  iNET  data  packet  structure  which  have  been  transmitted  over  an  aeronautical  channel.  This 


equalizer  is  especially  effective  when  the  known  data  bits  in  the  form  of  the  preamble  and  asynchronous 
marker  (ASM)  bits  contained  in  the  iNET-packet  structure  are  used  to  provide  an  alternative  method  of 
initializing  the  CMA  equalizer.  Since  this  decreases  the  convergence  time  of  the  adaptation  and  actually 
produces  a  lower  final  bit  error  count  when  compared  to  the  usual  method  of  center  tap  initialization. 
For  real-time  processing,  the  block  processing  CMA  algorithm  is  applied  to  a  data  stream  of  multiple 
packets,  where  the  bit  error  analysis  is  done  after  convergence  of  the  equalizer  is  achieved,  and  we  then 
continue  to  update  the  CMA  weight  vector  using  each  new  incoming  data  packet.  The  effect  of  the 
equalizing  CMA  radius  on  the  BER  performance  is  also  investigated. 


COMMUNICATION  SYSTEM  MODEL 

The  INET  packet  structure  is  made  up  of  the  known  preamble  of  length,  128  bits  and  ASM  is  of  length, 
64  bits,  together  with  the  actual  data  of  length,  6144  bits.  The  total  packet  length  is  6336  bits  and  this 
packet  structure  is  shown  in  Fig.  1.  The  transmitted  signal  is  SOQPSK-TG,  a  partial-response 
continuous  phase  modulation  (CPM)  with  a  constrained  ternary  alphabet.  The  precoder,  modulation 
index  and  frequency  pulse  are  described  in  [5,6]. 


Preamble 


ASM 


Data 


Fig.  1 :  INET  Packet  structure 


The  system  communication  model  is  shown  in  Fig.  2.  A  bit  stream,  formatted  as  shown  in  Fig.  1,  is 
modulated  using  SOQPSK-TG.  The  SOQPSK-TG  signal  experiences  multipath  propagation,  in  the  form 
of  an  LTI  system  with  impulse  response  h(n),  with  the  addition  of  white  Gaussian  noise.  The  received 
signal  is  equalized  by  an  FIR  filter  with  impulse  response  w(n).  The  equalizer  output  is  input  to  a 
symbol-by-symbol  SOQPSK-TG  demodulator  [6]  to  produce  the  bit  estimates  used  to  assess  the 
performance  of  the  equalizer.  Because  of  the  known  preamble  and  ASM  bits  in  the  transmitted  data 
packet,  it  is  possible  to  compute  the  minimum  mean  square  error  (MMSE)  equalizer,  which  is  then  used 
to  initialize  the  CMA  equalizer  for  efficiency.  The  CMA  equalizer  is  updated  by  using  each  received 
data  block  after  convergence  of  the  adaptation,  and  bit  error  rates  are  determined. 


v(n) 


Fig.  2:  SOQPSK  Communication  System  Model 


The  channel  is  a  single-input  single-output  (SISO)  system  characterized  by  the  input/output  relationship: 


x(n)  =  T,  h(k)s(n  -k)  +  v(n), 

k=-K , 


(1) 


where  h(n)  is  the  channel  impulse  response  which  is  non-causal  and  of  length  (K\  +K2  +1),  s(ri)  are 
samples  of  the  complex-valued  low-pass  equivalent  signal  corresponding  to  the  SOQPSK-TG 
modulated  carrier  within  the  packet  structure,  and  v(n)  is  AWGN  noise.  The  sequence  x(n)  comprises  N 
samples.  Assuming  the  received  signal  is  sampled  at  a  rate  equivalent  to  2  samples/bit,  N  is  at  least 
12672.  (The  actual  length  depends  on  the  length  of  the  channel.)  For  our  experiments  the  equalizer 
update  is  done  using  a  single  packet  of  transmitted  databits,  which  was  randomly  generated,  these  are 
modulated  and  transmitted  over  a  set  of  aeronautical  test  channels  [4].  The  equalizer  output  is  given  by 

y(n)=Ywn(l)x(n-l),  (2) 

/=-L, 


where  wn  =  \wn(-Ll) . w„(-l)  w„  (()) . wn(L2)\'  is  the  equalizer  weight  vector  at  iteration,  n  of  the 

CMA  adaptation  process. 


CMA  EQUALIZATION  ALGORITHM 

Because  SOQPSK-TG  is  constant  modulus,  it  is  expected  that  it  will  be  well-suited  for  equalization 
using  CMA.  The  MMSE  equalizer  is  the  filter  that  minimizes  the  mean  square  error  between  the 
equalizer  output  y(n)  and  the  channel  input  s(n)  [7].  By  initializing  CMA  with  the  MMSE  equalizer, 
there  will  be  no  phase  ambiguity  for  the  CMA  equalized  symbols. 

Description  of  CMA  Equalization  Algorithm 

The  constant  modulus  algorithm  (CMA)  is  a  blind  equalization  algorithm,  which  provides  an  adaptation 
based  on  a  specific  cost  function.  The  cost  function,  JCMA(y(n )) ,  which  is  minimized,  is  a  function  of  the 
distance  of  the  equalizer  output  from  a  circle  of  known  radius.  This  radius  is  detennined  from  the 
modulation  used  for  signal  transmission.  The  CMA  cost  function  is  given  by 

JCMA{y{n))  =  E^y{n)\2-R2)^,  (3) 

where  7G  is  the  radius  squared  of  the  CMA  circle  on  which  the  received  samples,  x(n)  lie.  This  cost 
function  performs  equalization  by  taking  into  account  the  distance  between  the  equalizer  output  samples 
and  the  pre-specified  radius  of  the  desired  signal.  For  SOQPSK-TG,  the  transmitted  signal  is  CPM  and 
of  the  form:  s(n)  =  exp  (j</>(n))  ,  so  the  value  of  CMA  radius  is  7G  =  1-  Note  that  because  SOQPSK-TG  is  a 
CPM  it  has  constant  modulus,  and  the  equalizer  operates  on  samples  of  the  complex-valued  low-pass 
equivalent  signal,  which  are  the  received  samples  from  a  transmitted  signal  that  was  constant  modulus. 
This  motivates  the  use  of  CMA. 


The  update  of  the  equalizer  weights  for  this  cost  function  is  based  on  a  stochastic  gradient  descent  rule 
given  by 


™n+ 1  =™„-M  VwJCMA(y(n ))  (4) 

where  p  is  the  algorithm  step-size,  the  equalizer  output  is  a  block  of  data  detennined  by  y(n)  =  wn'  xn , 
for  equalizer  weight  vector  wn ,  and  a  received  data  vector,  xn  =\x (n  +  . x(n ) . x(n-L2)]T  . 

The  gradient  vector  of  the  cost  function  is  given  as: 

V^CvOO)  =  e\*{  \y(n)f  -R2)  y(n)  x/j  •  (5) 

The  *  operator  denotes  the  complex  conjugate,  and  for  these  experiments  an  adaptive  step-size,  u  was 
used  in  (4)  which  is  given  (with  a=0.4)  by 

T 

X„  w/7 

**  ||x J  V 

In  general  the  CMA  equalizer  weight  vector  is  updated  starting  from  an  initial  vector,  wo  which  is 
center-tapped.  Previous  work  in  [3]  has  shown  the  improvement  in  performance  over  the  CMA  equalizer 
with  center-tap  initialization  of  the  CMA  equalizer  when  it  is  initialized  using  the  MMSE  equalizer. 
Since  the  MMSE-initialized  CMA  equalizer  does  converge  faster  to  achieve  a  lower  bit  error  rate.  For 
SOQPSK  the  value  of  the  CMA  radius  to  be  used  should  be  unity,  but  a  non-unity  radius  of  2  is  tested  in 
our  experiments  below  in  order  to  equalize  the  received  SOQPSK  modulated  signals  and  it  is  observed 
that  when  the  value  of  R2  in  equation  (3)  is  chosen  to  have  the  value  of  2,  instead  of  unity,  the  CMA 
equalization  works  much  better. 


Aeronautical  Telemetry  Channels 

Four  aeronautical  test  channels  with  impulse  responses  which  have  been  previously  measured  from 
channel  sounding  experiments  at  Edwards  AFB  [4]  were  used  in  our  experiments.  Characteristics  of 
these  four  channels  used  are  provided  in  Table  1  in  terms  of  total  channel  delay  spread  and  total  number 
of  non-zero  multi-paths.  These  channels  represent  snapshots  of  the  flight  path  of  an  aircraft.  Channels  1 
and  2  were  captured  on  a  taxiway  ,  channel  4  was  captured  during  takeoff  on  runway  22  L,  and  channel 
6  is  airborne  flight. 


Table  1 :  Channel  Information 


Channel  # 

1 

2 

4 

6 

Channel  Length 

9 

20 

19 

4 

No.  of  Non-zero  taps 

3 

8 

9 

4 

EXPERIMENTAL  SETUP 


Implementation  of  the  CMA  for  PAQ 

For  this  implementation  the  CMA  equalizer  weight  vector,  wo  is  initialized  using  the  minimum  mean 
square  (MMSE)  equalizer,  which  is  computed  using  the  known  (the  preamble  and  ASM)  databits 
contained  in  the  data  packet.  For  this  application  CMA  is  to  be  evaluated  on  a  real-time  system  which 
consists  of  data  stream  received  as  multiple  packets.  Thus  for  the  case  of  real-time  processing,  there  is 
an  initialization  phase  during  which  the  incoming  data  packets  are  being  processed  to  achieve 
convergence  of  the  CMA  equalizer.  This  phase  may  take  50  to  100  data  packets,  after  which  the  bit  error 
analysis  can  be  done  to  count  total  errors  on  the  following  data  packets.  Two  real-time  approaches  are 
considered  here  as  described  in  the  simulation  procedures  below. 

Simulation  Procedure  1  (Multiple  Packet  -  Adaptive  startup) 

This  is  referred  to  as  the  Block-by-Block  (BxB)  update  procedure. 

Step  1:  Generate  an  iNET  packet  at  2  samples/bit.  Send  the  packet  through  the  test  channel  and  add 
noise  to  create  the  received  packet. 

Step  2\  Initialize  the  equalizer  filter. 

Step  3:  Equalize,  update  the  equalizer  filter  coefficients  based  on  the  received  packet  produced  in  Step  1. 
Go  to  step  1  until  the  equalizer  has  converged  (50  packets),  then  proceed  to  Step  4. 

Step  4\  Generate  a  new  iNET  packet  at  2  samples/bit.  Send  the  packet  through  the  test  channel  and  add 
noise  to  create  a  received  packet. 

Step  5:  Equalize,  update  the  equalizer  filter  coefficients,  and  count  the  number  of  bit  errors. 

Step  6:  Repeat  Step  4  and  Step  5,  500  times. 

Simulation  Procedure  2  (Single  Packet  -  Static  startup) 

Step  1:  Generate  an  iNET  packet  at  2  samples/bit.  Send  the  packet  through  the  test  channel  and  add 
noise  to  create  the  received  packet. 

Step  2:  Initialize  the  equalizer  filter. 

Step  3:  Iterate  on  the  received  packet  produced  in  Step  1,  updating  the  equalizer  filter  once  per  packet, 
until  the  equalizer  has  converged  (50  iterations). 

Step  4\  Generate  a  new  iNET  packet  at  2  samples/bit.  Send  the  packet  through  the  test  channel  and  add 
noise  to  create  a  received  packet. 

Step  5:  Equalize,  update  the  equalizer  filter  coefficients,  and  count  the  number  of  bit  errors. 

Step  6:  Repeat  Step  4  and  Step  5,  550  times. 

The  main  difference  between  these  two  procedures  is  as  follows:  For  Simulation  Procedure  1,  the  filter 
update  equation  (4)  achieves  convergence  by  updating  from  startup  using  multiple  received  data  packets. 
For  Simulation  Procedure  2  at  startup,  a  single  fixed  data  packet  is  used  for  the  update  of  equation  (4) 
until  convergence  at  50  iterations,  and  only  after  this  does  its  update  continue  with  a  new  data  packet 
being  processed  at  each  iteration. 

In  the  section  below  we  present  results  for  the  two  real-time  processing  schemes  outlined  above, 
and  we  compare  performance  of  this  MMSE-initialized  CMA  equalizer  for  the  communication  system 
model  which  was  described  in  Section  II  based  on  the  test  channels  outlined  in  Table  I.  We  consider  the 


effect  of  changing  some  parameters  in  order  to  achieve  the  best  performance  indicated  by  the  baseline 
bit  error  rates  from  [3],  and  in  particular  we  consider  the  effect  of  using  a  CMA  radius  squared  with  a 
value  of  2  for  SOQPSK  equalization.  Note  that  the  baseline  performance  is  determined  by  running  steps 
1-3  of  procedure  2,  and  averaging  over  a  500  Monte  Carlo  simulations.  In  addition  we  also  compare  the 
effect  of  using  a  different  CMA  radius  squared  value. 

PERFORMANCE  EVALUATION  RESULTS 

In  this  section  we  present  results  which  show  the  performance  of  the  MMSE-initialized  CMA  equalizer 
for  the  communication  system  model  described  above  using  four  of  the  aeronautical  telemetry  channels. 
Since  the  CMA  equalizer  is  initialized  with  the  MMSE  equalizer,  the  MMSE  equalizer  performance  is 
included  in  all  plots  as  well  as  the  baseline  performance  from  [3],  and  the  theoretical  lower  bound  curve 
for  AWGN  channel.  We  compare  the  bit  error  curves  of  the  multiple  packet  startup  procedure  1  to  the 
single  packet  startup  procedure  2,  where  the  initialization  for  both  procedures  is  done  over  50  data 
packets,  and  the  bit  error  analysis  is  done  by  counting  errors  over  600  packets  of  received  data. 


Fig.  3:  Single  Packet  Processing  vs.  Block-by-Block  CMA  Equalizer  for  Channel  1 

R2  =  2  (left)  versus  Rj  =  1  (right) 


Fig.  4:  Single  Packet  Processing  and  Block-by-Block  CMA  Equalizer  for  channel  3 

R2  ~  2  (left)  versus  Rj  =  1  (right) 


Channel  4 


Fig.  5:  Single  Packet  Processing  and  Block-by-Block  CMA  Equalizer  for  channel  4 

R2  ~  2  (left)  versus  Ri  =  1  (right) 


Fig.  6:  Single  Packet  Processing  and  Block-by-Block  CMA  Equalizer  for  channel  6 

R-2  =  2  (left)  versus  R/  =  1  (right) 


In  Figures  3  -  6  we  compare  the  performance  of  the  CMA  equalizer  using  the  multiple  packet  startup 
(Simulation  Procedure  1)  to  the  single  packet  startup  (Simulation  Procedure  2)  for  two  different  values 
of  the  CMA  radius  squared  fC.  Radius  squared  R2  =  1  is  on  the  right,  while  R2  =  2  is  the  plot  on  the  left 
of  each  numbered  figure.  We  note  that  for  both  of  the  procedures  tested,  the  bit  error  rate  perfonnance  is 
generally  better  for  R2  =  2  than  for  /C  =  1 ,  especially  at  the  higher  SNRs.  In  fact  at  low  SNR,  the 
perfonnance  of  the  MMSE-initialized  CMA  using  radius  R2  =  1  is  worse  than  the  MMSE  equalizer, 
while  the  performance  of  these  two  procedures  is  better  than  MMSE  when  the  CMA  radius  squared  R2  = 
2  is  used.  Also  from  the  curves  on  the  right  (R2  =1)  neither  of  these  two  procedures  appears  to  have  an 
bit  error  rate  advantage  over  the  other.  It  has  been  observed  that  as  the  total  number  of  packets  processed 
increases,  the  error  rate  for  the  single  packet  startup  (Simulation  Procedure  2)  tends  to  increase  more 
rapidly  than  that  of  multiple  packet  startup  (Simulation  Procedure  1). 


CONCLUSIONS 


In  this  work,  we  have  presented  two  block-processing  CMA  equalization  schemes  for  SOQPSK-TG 
modulated  data,  which  has  been  transmitted  over  an  aeronautical  channel.  Since  the  transmitted  data 
packets  contain  known  data  for  the  telemetry  application  of  interest  (iNET  data  packet  structure)  the 
MMSE  equalizer  has  been  used  as  an  initializer  for  CMA  in  this  application.  We  have  implemented, 
tested  and  compared  the  use  of  a  multiple  packet  startup  (Simulation  Procedure  1)  and  a  single  packet 
startup  (Simulation  Procedure  2)  for  the  CMA  algorithm.  In  general,  the  MMSE-initialized  CMA 
equalizer  using  radius  squared  Ri  =  2  achieves  better  performance  than  using  radius  1  for  both 
procedures  introduced  here.  In  future  work  we  will  focus  on  real-time  implementation  issues,  to  consider 
whether  there  are  conditions  under  which  a  re-initialization  of  the  equalizer  may  be  required  using  the 
known  databits  contained  in  each  packet.  The  perfonnance  of  the  multiple  packet  startup  procedure  1  is 
generally  more  stable  producing  a  lower  error  rate  than  single  packet  startup  case.  It  is  unclear  exactly 
why  the  non-obvious  choice  of  CMA  radius  2  provides  lower  bit  error  rates,  and  this  continues  to  be  a 
topic  of  further  investigation. 
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